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METHODS AND SYSTEMS FOR GENERATING A TERRAIN ELEVATION MAP 
10 IN A CARTESIAN FORMAT 
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Priority Claim 

This application claims the benefit of U.S. Provisional Application Ser. No. 
60/528,297, filed December 10, 2003, which is hereby incorporated by reference. 

Field of the Invention 
20 This invention relates generally to navigation and, more specifically, to displaying 

terrain elevation maps in a vehicle. 

Background of the Invention 

Terrain elevation maps are frequently used in aviation for purposes of maintaining 
pilot spatial situation awareness, and for other purposes that include airborne radar signal 
25 processing. For display, or for processing purposes, the map of terrain is represented on a 
two dimensional planar display in an x/y Cartesian format, even though the terrain is 
superimposed on the spherically curved earth surface. Because the radius of curvature is 
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large, and the size of the region represented by the map is small with respect to the radius of 
the earth, the "flat earth" representation of the terrain elevation map is acceptable. 

The use of terrain elevation data in a Cartesian format is complicated by the fact that 
terrain data is typically provided as points uniformly sampled in latitude and longitude. It is 
5 often the practice to present terrain data in the latitude/longitude format as an approximation 
to the desired x/y format. This representation of terrain elevation in the latitude/longitude 
format is reasonably good at latitudes near the equator. However, terrain elevation data 
represented in the lat/long format presents some processing anomalies at other latitudes. 

Typically, a terrain elevation map in an aircraft, or in any moving vehicle, represents 

10 the terrain over a limited size region in the vicinity of the vehicle. The size of the region 
depends on display, or other, requirements. As the vehicle moves, it is necessary to update 
the map data to incorporate the new terrain. If the map is maintained with an orientation 
referenced to true north, the reference frame of a generated map has to rotate to maintain the 
orientation to true north as the vehicle travels. This will occur at any latitude except at the 

15 equator. As a result, if the map is maintained referenced to north, a map update that includes 
this reference frame rotation must consist of re-building the entire map in order to maintain 
the north reference. This is necessary despite the fact that new terrain entering the mapped 
region may constitute a relatively small portion of the map. Re-building the entire map can 
result in substantial computation burden and cost. 

20 Therefore, there exists a need to provide a terrain elevation map that does not 

introduce excessive distortion, and which permits an incremental update to the map as the 
vehicle moves that prevents the need to periodically re-build the map. 
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Summary of the Invention 

The present invention provides methods, systems, and computer program products for 
re-sampling terrain data in a vehicle from latitude/longitude format into a Cartesian format, 
to permit incremental updates of the terrain elevation map as the vehicle travels. 
5 In one embodiment the system includes memory for storing terrain elevation data, a 

navigation component for generating aircraft position information, and a processor in 
communication with the memory and the navigation component. The processor retrieves 
terrain elevation data from the memory based on the generated aircraft position information, 
and generates a Cartesian coordinate-referenced terrain elevation map based on the retrieved 
1 0 terrain elevation data and the aircraft position information. 

The system also includes a display processor in communication with the processor 
that generates a map image. 

In one aspect of the invention, the navigation component includes a ground proximity 
warning system that generates one or more of a ground proximity caution or warning based 
15 on the generated terrain elevation map. 

In another aspect of the invention, the system includes a radar system that generates a 
volumetric buffer based on the terrain elevation map and radar signals generated by the radar 
system. 

Brief Description of the Drawings 

20 The preferred and alternative embodiments of the present invention are described in 

detail below with reference to the following drawings. 

FIGURE 1 is a block diagram of an exemplary system formed in accordance with the 
present invention; 

FIGURES 2-4, 6, and 7 are flow diagrams of an exemplary process performed by the 
25 system shown in FIGURE 1 ; and 

FIGURE 5 is an image of reference frames on the earth. 
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Detailed Description of the Preferred Embodiment 

FIGURE 1 illustrates a block diagram of a vehicle 20, such as an aircraft, formed in 
accordance with an embodiment of the present invention. The vehicle 20 includes a 
processor 30 and a database 34 in data communication with the processor 30. The vehicle 20 
5 may include one or more of a navigation system 32, a display 36 having a display processor, 
and a radar system 38. The navigation system 32, the display 36, and the radar system 38 are 
all in data communication with the processor 30 or the database 34. 

The processor 30 generates a terrain elevation map by re-sampling terrain data stored 
in the database 34 from a latitude/longitude format into a Cartesian format based on aircraft 
10 position information received from the navigation system 32. The processor 30 also 
incremental updates the terrain elevation map as the vehicle 20 travels. The terrain elevation 
map is stored in the database 34 or transmitted to the navigation system 32, the display 36, or 
the radar system 38. 

In one embodiment, the processor 30 is a general-purpose digital computer that 
15 interacts with the navigation system 32, the database 34, the display 36, and/or the radar 
system 38 and operates according to processes described below. 

In one embodiment, the terrain elevation map is presented on the display 36. In 
another embodiment, the navigation system 32 uses the terrain elevation map for use by a 
ground proximity warning component or a navigation component, or for other associated 
20 purposes. In still another embodiment, the radar system 38 uses the terrain elevation map in 
various operations, such as radar tilt control or in generation of a volumetric radar return 
buffer. 

FIGURE 2 illustrates a flow diagram of an exemplary process 70 performed by the 
processor 30 (FIGURE 1). First, at block 74, the processor 30 retrieves terrain elevation data 
25 from the database 34 based on aircraft position information received from the navigation 
system 32. Next, at block 76, the processor 30 generates a Cartesian coordinate-referenced 
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terrain elevation map based on the retrieved terrain elevation data and aircraft position 
information received from the navigation system 32. Generation of the Cartesian coordinate- 
referenced terrain elevation map is described in more detail below in FIGURE 3. The terrain 
elevation map is stored in the database 34, block 78, sent to the display 36 with display 
5 processor for presentation, block 80, sent to the navigation system 32 for use in warning and 
navigation functions, block 82, and/or sent to the radar system 38 for use in radar signal 
processing, block 84. 

In an alternate embodiment, the terrain elevation map is transmitted from the 
database 34 to the navigation system 32, the display 36, or the radar system 38. 

10 As shown in FIGURE 3, a process 100 is shown for generating and updating the 

terrain elevation map (block 76, FIGURE 2). The process 100 begins at a block 104 where an 
initial Cartesian coordinate-referenced terrain elevation map is generated. Generation of the 
initial terrain elevation map is described in more detail below in FIGURE 4. Next, at a 
block 106, the Cartesian coordinate-referenced terrain elevation map is updated. Map 

15 updating is described in more detail below with regards to FIGURE 6. At a block 108, the 
processor 30 updates a rotation angle. Updating of the rotation angle is described in more 
detail below with regards to FIGURE 7. If at decision block 112 the processor 30 is still 
active (i.e. the processor 30 has be deactivated or requested to stop generating a terrain 
elevation map), the process 100 ends, otherwise, the process 100 returns to block 106. It can 

20 be appreciated that the process 100 may be terminated at any point. 

As shown in FIGURE 4, a process 130 for generating an initial terrain elevation map 
is illustrated. First, at a block 132, a transform matrix for transforming points in a terrain 
elevation map from a Earth-Centered Earth-Fixed (ECEF) frame to a Stabilized Local-Level 
(SLL) frame is set equal to a transform matrix for transforming points in a terrain elevation 

25 map from a Earth-Centered Earth -Fixed (ECEF) frame to a Local-Level North Referenced 
(LLNR) frame. FIGURE 5 illustrates a geometric view of an ECEF coordinates (x ece f 5 y e cef> 
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Zecef) and an LLNR coordinates (xn nr , yiinr, zn nr ). Next, at a block 136, a rotation angle 
between the LLNR frame and the SLL frame is set to zero. Next, at a block 138, a position 
(x e cef> Yecef, z ece f) is calculated for each point in the terrain elevation map in ECEF frame 
based on position information generated by the navigation system 32. At block 140, latitude 
5 and longitude (lat/long) are calculated for each point in the terrain elevation map in ECEF 
frame based on the calculated positions (x ece f, y eC ef, z ece f). Next, at a block 144, terrain 
elevation data is mapped into the terrain elevation map in SLL frame based on the calculated 
lat/long, the SLL transform matrix, and the stored terrain elevation data. 

As shown in FIGURE 6, a process 180 (block 106, FIGURE 3) for updating the 

10 terrain elevation map produced by the process 130 is shown. First, at a block 184, 
incremental distance measurements in x and y direction of the SLL frame are determined 
based on information from the navigation system 32. At decision block 186, the processor 30 
determines if the vehicle has moved an incremental x or y distance more than a threshold 
amount. If the vehicle has not moved an incremental x or y distance more than a threshold 

15 amount, then the process 180 returns to block 184. If the vehicle has moved an incremental x 
or y distance more than a threshold amount, a corresponding row or column is added to the 
terrain elevation map based on the present SLL transform matrix and the new row or column 
is populated with terrain elevation data. 

The terrain elevation map is maintained in the SLL frame in a Cartesian (i.e., x/y) 

20 format. In one embodiment the terrain elevation map is maintained with the vehicle (e.g., 
aircraft) near the center. If the vehicle moves a distance equal to one cell in x or y (or actually 
in the curvilinear coordinates xi and eta), then a row or column of new terrain data is needed 
and a row or column is dropped off of the terrain elevation map behind the vehicle based on 
the direction of travel. When adding a new cell, the curvilinear coordinates of the points 

25 (cells) in the new row or column are determined. The curvilinear coordinates represent the 
distance along the surface of the earth from the reference point (i.e., where the SLL z-axis 
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penetrates the surface of the earth). The curvilinear coordinates are converted into rotation 
angles by dividing by the radius of the earth. The position of the point in the SLL frame is 
calculated using the rotation angles. The position is transformed to the ECEF frame by 
multiplying by the SLL-to-ECEF transformation matrix. The transformed position is used to 
5 calculate the lat/long of the point. The terrain elevation data associated with the calculated 
lat/long is then entered into, assigned to, or associated with the point. 

As shown in FIGURE 7, a process 200 is shown for updating the rotation angle about 
the vertical direction (block 108, FIGURE 3). First, at a block 204, vehicle incremental 
displacement values are determined. At a block 206, the SLL/ECEF transform matrix is 

10 updated based on the determined incremental displacement values. Then, at block 210, the 
rotation angle between SLL and LLNR frames is updated. 

The rotation angle is used for performing geometric calculations with the terrain 
elevation map, such as those done for the Volumetric Buffer radar processing and for display. 
For display, the terrain elevation map is preferably oriented with respect to true north. The 

15 rotation angle indicates how much to rotate the terrain elevation map when processing it for 
display. 

Mathematical Description 

Locations in a map plane include two-dimensional coordinates (£ 77), where £ and t] 
represent distances in a display plane in orthogonal directions. On a curved earth surface, the 
20 (£ 77) coordinates correspond to arc lengths along the mean surface of the earth. Therefore, 
the (£ 77) coordinates map into angle coordinates {</>, 0) via the following relations: 

0=^-, (i) 

Y Y 

earth earth 

25 where Y earth is the mean earth radius. 
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A z-axis points in the local vertical direction, a position vector [x,y,z] T of a point in 
the display plane defined by the coordinates (£ rj) is calculated as follows: 



X 




cos 6 sin ^ 


y 




sin# 


z 




cos 0 cos <f> 



(2) 
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where (j> and 0are obtained from £ and rj using Equation (1). The position defined in 
Equation (2) is normalized by the radius of the earth. From Equation (2) £ signifies an arc 
length in the x direction, and rj signifies an arc length in the y direction. 

A reference frame where the z-axis is in the local vertical direction is selected. An 
example of a commonly used reference frame is the local-level north referenced (LLNR) 
frame. FIGURE 5 illustrates the orientation of the LLNR with respect to an earth-centered 
earth-fixed (ECEF) reference frame. 

A position vector in the ECEF frame is obtained from the following. 





X 




cos 0 sin O 


15 


y 




sin© 




z 


ECEF 


cos© cos O 



(3) 



where 0 is the latitude, and O is the longitude. The latitude and longitude are 
calculated as follows: 



20 



0 = sin y ECEF , O = arctan(z £C£F , x ECEF ) 



(4) 



where x E cef, yECEF, and z E cef are the components of the position vector in the ECEF 
frame, and arctan(x,y) is a four quadrant arctangent function. 
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At this point, we can calculate position vectors in the LLNR frame from curvalinear 
coordinates (£ rj) using Equations (1) and (2), and we can calculate latitude and longitude 
from the position vector defined in the ECEF frame. What is missing is a means of 
transforming the position vector defined in the LLNR to a position vector defined in the 
ECEF. We can accomplish this by performing the following operation: 



x 

y 

L ^ J ECEF 



= A 



LLNR I EF 



(5) 



LLNR 



where [x, y, z\ E cef describes the location of a point in the ECEF frame, [x, y, z\ LLN r 



10 describes the location in the LLNR frame and A L lnr/ef is the transformation matrix given by 

cos<D - sin© sin O cos© sin O 

0 cos© sin© (6) 

-sin d> - sin© cos O cos© cos O 



LLNR I EF 



where 0 is the latitude of the vehicle, and <D is the longitude of the vehicle. 

Databases of terrain elevation are provided in the form of sample points with respect 
15 to a reference frame fixed to the earth. These samples are provided as samples with some 
spacing in latitude and longitude. In order to get a terrain elevation map with a proper spatial 
relationship between sample points, we need to be able to calculate latitude and longitude for 
points sampled in (£ 77) space. 

The operations to perform this task are as follows. For each of the points sampled in 
20 (£ rj) , the application of Equations (1) and (2) results in the position vector in the LLNR. 
The position vector is then multiplied by the transformation matrix (Equations (5) and (6)) to 
give the position vector in the ECEF frame. These coordinates are then applied to Equation 
(4) to get the latitude and longitude of the sample point. Using latitude and longitude, the 
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terrain elevation is extracted from the terrain database. The result is a terrain elevation map 
with the proper spatial relationships between sample points. 

In an alternate embodiment, the processor 30 augments the terrain elevation map only 
along its leading edges as the vehicle moves, thereby, eliminating regenerating the entire 
5 terrain elevation map. This embodiment ensures that the rows and columns of the map that 
are generated are valid (i.e., sufficiently accurate) during the entire time it takes for the 
vehicle to traverse the map. However, performing this type of incremental update using the 
LLNR reference frame is not practical since maintaining the orientation with true north 
requires the rotation of the reference frame, and so the rotation of the rows and columns of 

10 the map. However, maintaining the terrain map in a different reference frame can allow the 
desired incremental updating. 

This new reference frame is similar to the LLNR in that the z-axis is in the local 
vertical direction. However, the new reference frame is not referenced to north, but instead is 
in a sense "stabilized" so that there is no rotation. This new reference frame will be referred 

15 to as the Stabilized Local Level (SLL) frame. 

The use of the SLL frame is accomplished by maintaining the SLL-to-ECEF 
transformation matrix. Then the computation of terrain elevation map points is performed in 
a manner similar to the procedure described above, but with the SLL-to-ECEF transform 
matrix used in place of the LLNR-to-ECEF transform. 

20 The SLL-to-ECEF transform matrix A S luef is initialized to be equal to A L lnr/ef, 

which is calculated from the vehicle's latitude and longitude using Equation (6). However, 
whereas the A LLNR/EF matrix is updated at some time interval by recalculation using the 
vehicle's new latitude and longitude, the update of the SLL-to-ECEF transform is 
accomplished by performing an incremental rotation of the SLL frame by rotation angles 

25 defined in the SLL frame. These rotation angles are obtained by accumulating the vehicle 
distance traveled since the last update using data from the navigation system 32. This results 
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10 



in an integral displacement of 5% and Stj defined in the SLL frame. 8% and 8tj are easily 
calculated using data from a navigation system. However, because the SLL is incrementally 
rotated with respect to the LLNR, the calculation of <5£ and 8rj from navigation system data 
must take into account the rotation angle. Initially, the incremental rotation angle is zero, 
because the SLL is initially the same as the LLNR. As the vehicle travels, the rotation angle 
deviates from zero, so it is necessary to calculate the rotation angle as part of the processing. 
The calculation of the rotation angle is described later in this document. 

Given the incremental displacements 8% and Srj since the last update, the incremental 
rotation angles for the update are calculated by dividing by the radius of the earth, giving: 

(7) 



'earth 



'earth 



These incremental rotations are then used to rotate the SLL frame,- resulting in a new 
SLL-to-ECEF transform obtained from the following calculation. 



15 



{A SLL I EF )next ~ ^AsLL I EF ) 



previous 



cos 8$ -sin 86 sin 8$ cos 89 sin 8</> 

0 cos£0 sm89 

-sin<S# -sin 89 cos 8<f> cos 89 cos 8<j> 



(8) 



While the LLNR frame is forced to have the z-axis directed through the vehicle 
position in the local up direction, it is possible for errors to build up as the process given in 
Equation (8) is iteratively applied. The build-up of errors results if the SLL z-axis drifting 
20 away from the vehicle position. However, the SLL frame can be "locked" to the vehicle by 
calculating the following difference vector after each update: 





' 0 " 




" 0 ' 


A - A T A 

SLL/ EF LLNR 1 EF 


0 




0 




J earth _ 




J earth _ 
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The matrix product A 1 sll/efAllnr/ef gives the LLNR-to-SLL transform. Therefore, the 
vector A is equal to the vector difference between the z-axis in the LLNR (which points to 
the vehicle position) transformed to the SLL, and the z-axis in the SLL. This is essentially the 
5 residual vehicle position error in the SLL after the update. The SLL is kept "locked" to the 
vehicle by setting the aircraft incremental displacement after the update to values derived 
from A as follows: 8E, = the x component of A, 8r\ = the y component of A. Then the 
navigation system data is used to accumulate additional contributions to these incremental 
displacements until the next SLL frame update. 

10 The navigation system 32 typically provide heading and velocity vector data 

referenced to true north. However for use with the SLL, heading and track angles referenced 
to a principal axis in the SLL frame (the y-axis for example) are used. Therefore, the rotation 
angle between the SLL and LLNR frames are determined for transforming velocity and 
heading in the LLNR (from navigation sensors) to velocity and heading in the SLL. In one 

15 embodiment, the rotation angle between the SLL and LLNR frames is obtained by 
calculating the orientation of north (i.e., the LLNR frame y-axis) in the SLL frame as 
follows: 



20 



= A T A 

SLL I EF LLNR I EF 



SLL 



(10) 



The rotation angle is obtained from the x and y components by 

r=**ctm(x su> y su ) (11) 
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Then to determine the corresponding angle to use for navigation within the SLL 
frame, the angle y should be added to the north-referenced angle obtained from the 
navigation system, so that 



5 



<f>SU =<t>LLNR +7 



(12) 



where feu is a bearing angle (i.e., heading or track angle) in the SLL frame and fe Lm 
is the corresponding true north referenced bearing obtained from a navigation system. The 
rotation angle ^is calculated at the update time and used during the period between updates. 



as noted above, many changes can be made without departing from the spirit and scope of the 
invention. For example, the steps performed in the processes above may be performed in 
various order without departing from the scope of the invention. Accordingly, the scope of 
the invention is not limited by the disclosure of the preferred embodiment. Instead, the 
1 5 invention should be determined entirely by reference to the claims that follow. 
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While the preferred embodiment of the invention has been illustrated and described, 
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